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Abstract The mammalian spatial navigation system makes use of several different sensory information 
channels. This information is then converted into a neural code that represents the animal’s current 
position in space by engaging place cell, grid cell, and head direction cell networks. In particular, sensory 
landmark (allothetic) cues can be utilized in concert with an animal’s knowledge of its own velocity (id- 
iothetic) cues to generate a more accurate representation of position than (idiothetic) path integration 
provides on its own (Battaglia et al, 2004). We develop a computational model that merges path integra¬ 
tion with information from external sensory cues that provide a reliable representation of spatial position 
along an annular track. Starting with a continuous bump attractor model, we allow for the possibility of 
synaptic spatial heterogeneity that would break the translation symmetry of space. We use asymptotic 
analysis to reduce the bump attractor model to a single scalar equation whose potential represents the 
impact of heterogeneity. Such heterogeneity causes errors to build up when the network performs path 
integration, but these errors can be corrected by an external control signal representing the effects of 
sensory cues. We demonstrate that there is an optimal strength and decay rate of the control signal 
when cues are placed both periodically and randomly. A similar analysis is performed when errors in 
path integration arise from dynamic noise fluctuations. Again, there is an optimal strength and decay of 
discrete control that minimizes the path integration error. 

Keywords neural field, sensory feedback, spatial navigation, stochastic differential equation 


1 Introduction 

Animals have a remarkable ability to accurately navigate over large distances (Geva-Sagiv et al, 2015). 
For instance, birds can utilize neural systems that sense the earth’s magnetic field, orienting themselves 
geocentrically (Cochran et al, 2004; Wu and Dickman, 2012). This is in contrast to the systems studied in 
mammalian species, who are typically shown to use path integration (Etienne et al, 1996; McNaughton 
et al, 2006). Path integration models of spatial navigation assume mammals have knowledge of their 
direction and speed of motion, which networks of the brain can then integrate to encode the path of 
their idiothetic motion (Samsonovich and McNaughton, 1997). However, there are a number of potential 
sources of error to this mechanism. The nervous system itself is prone to a wide variety of noise sources due 
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to channel fluctuations, synaptic failures, or even stochastic network-wide events (Faisal et al, 2008). This 
could lead to faulty communication of velocity or head direction signals, or it could corrupt the storage 
of the position signal. Furthermore, the network that integrates the velocity signal may be comprised of 
architecture that is heterogeneous, providing an imperfect summation of velocity inputs (Brody et al, 
2003). Any inaccuracy in the represented position or velocity will be compounded over time, as the 
neural code continues to trace the animal’s true path (Knierim et al, 1995; Valerio and Taube, 2012). 
Fortunately, path integration is not the sole navigational technique of the mammalian brain; landmarks 
detected by the sensory system help anchor and correct the integrated velocity signal (Fig. lA) (Collett 
and Graham, 2004; Etienne et al, 1996; Solstad et al, 2008). 

Several experiments have demonstrated that mammals’ representation of space is sharpened in the 
presence of sensory cues (Aikath et al, 2014; Battaglia et al, 2004; Ulanovsky and Moss, 2011; Zhang 
et al, 2014). Experiments typically compare place fields of individuals cells - spatial locations where the 
cell becomes active - in the presence and absence of sensory landmarks (e.g., steel brush or ticking clock; 
Eig. IB). Eor instance, Battaglia et al (2004) recorded from hippocampal place cells in rats moving on 
an annular track. When there were no sensory cues along the track, the place fields of individual cells 
differed substantially, depending on whether the rat was moving clockwise or counterclockwise around the 
annulus. This suggests there was some drift in animals’ neural representation of their position. However, 
when several position landmarking cues were placed along the track, the clockwise and counterclockwise 
place fields of individual cells were considerably correlated. This suggests the sensory cues tightened the 
navigation systems fine representation of the animal’s spatial position. Similar effects have been observed 
in brown bats, whose echolocation signals provide a brief burst of rich sensory information, sharpening 
the animal’s place fields (Ulanovsky and Moss, 2011). Thus, sensory cues appear to provide a correction 
mechanism for the many sources of error that disrupt position representation and broaden place fields. 

We bring these two features of the mammalian navigation system - path integration and sensory 
feedback - together into a single model. Our main focus is the role sensory feedback can play in cor¬ 
recting the path integration signal. Errors in the path integration signal will arise, in our model, due to 
internal disruptions of an accurately delivered velocity input. The model represents position by utilizing 
a perturbation of a continuous bump attractor (Amari, 1977; Zhang, 1996). Bump states arise in these 
models due to a combination of strong local excitation and broadly tuned inhibitory feedback (Wang, 
1999). In translation symmetric networks, bumps can be formed with their center of mass at any location. 
However, such well-balanced architecture is unlikely to occur in actual networks of the brain, which tend 
to be spatially heterogeneous (Brody et al, 2003; Renart et al, 2003). Such symmetry-breaking in the 
bump attractor network leads to system states that drift toward a finite number of discrete attractors, 
so the long term dynamics are weakly correlated with the input signal (Itskov et al, 2011; Kilpatrick 
et al, 2013; Zhang, 1996). Eurthermore, any dynamic fluctuations in the voltage or synaptic signals of 
the network can lead to diffusive wandering of the bump state that will also degrade its signal represen¬ 
tation (Burak and Eiete, 2012; Compte et al, 2000; Kilpatrick and Ermentrout, 2013). Our study mainly 
focuses on how external control can reduce the deleterious effects of both spatial heterogeneity and noise, 
ultimately improving long term accuracy of the network’s path integration. 

Our model is intended to describe the neural activity of place cell networks in the hippocampus 
(O’Keefe and Burgess, 1996). Based on a simplified version of the model by Samsonovich and McNaughton 
(1997), neural activity bumps are propelled around the network by external velocity inputs that introduce 
spatial asymmetry into the balance of excitation and inhibition. A similar mechanism was utilized by 
Burak and Eiete (2009) in a two-dimensional model of grid cell activity. Thus, this framework is a well 
accepted model of position encoding cells in hippocampus, entorhinal cortex, and the vestibular system 
(Zhang, 1996). Since we are modeling motion along an annulus (Eig. IB), we have restricted the network 
to a one-dimensional periodic domain. Sensory cues are assumed to provide a reliable estimate of the 
animal’s true current position. This position is then compared with the place cell network’s estimate of 
position. Any discrepancy in the position estimate is then translated into a corrective velocity input, 
which is added to the baseline velocity input (Eig. lA). Even when the cues occur discretely in space, 
this mechanism works well for reducing the long term error in the position estimate. 

In section 2, we introduce the neural field model of spatial navigation, which combines path integration 
and sensory cue feedback. Next, we derive a low-dimensional approximation for the dynamics of bump 
position in the neural field model (section 3). This reduction reveals the relative influence of velocity 
inputs, sensory feedback, and heterogeneity on the animal’s perceived position of its current location. 
Ultimately, this allows us to calculate the impact of various control strategies on the error between the 
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Fig. 1 Mammalian spatial navigation network with sensory feedback. The animal utilizes its own velocity to update its 
remembered position (path integration) and corrects this memory with sensory cues that serve as position landmarks. 
(A) Schematic of the underlying neuronal network demonstrates the place cell network receives direction input from the 
animal’s velocity signal v(t). Its position estimate is compared with the true position read-out from a sensory cue, and 
this error is then used to generate a control input signal Vc{t) back into the place cell network. (B) Illustration of the 
experiments by Battaglia et al (2004), showing an annular track with various object placed in the environment to provide 
the animal with sensory cues. 


animars perceived position and true position (section 4). Our main finding is that there is an optimal 
control strength at which the long term error of the network is minimized. Our findings were similar 
in the case that errors arose due to dynamic noise fluctuations (subsection 4.2), rather than synaptic 
heterogeneities (subsection 4.1). In this case, the low-dimensional approximation of the neural field is a 
stochastic differential equation whose variance we can evaluate explicitly. 


2 Sensory control in velocity-integrating place cell networks 


We employ a neural field model of velocity integration that sustains a bump attractor of neural activity 
in the absence of any inputs. Amari (1977) pioneered the scalar neural field model as a reduction of the 
excitatory-inhibitory model of Wilson and Cowan (1973), but the incorporation of velocity inputs that 
shift the bump around the spatial domain is more recent. Originally developed as a model of the head 
direction system (Zhang, 1996), velocity-integrating networks introduce an external input that alters the 
shape of the recurrent architecture (McNaughton et al, 1991). As a result, a moving bump, rather than a 
stationary bump, becomes the stable solution to the model equations. This model has since been extended 
to account for place fields and grid cell fields in planar systems (Burak and Fiete, 2009; Samsonovich 
and McNaughton, 1997). The fully general form of our neural field model is given 


du{x, t) = — 


f 

J —71 


u{x,t)+ w{x,y)f{u{y,t))dy+ v 


it) r 

J —7T 


Wv{x - y)f{u{y,t))dy 


dtedW{x, t) (2.1) 


where u{x^ t) denotes the total synaptic activity at a position x G [—tt, tt] at some point in time t. While 
X labels the position of neurons in the network it also corresponds to location in the environment, so the 
domain Q = [—7r,7r] is taken to be periodic as it represents an annular track (Fig. IB). 

The function w{x,y) represents the synaptic connectivity between neurons, which we model as a 
translationally symmetric unimodal function icq, modified by spatial heterogeneity or odd asymmetry 
0, so 

w{x, y) := (1 + awu{y))wo{x - y - 4>) (2.2) 

with cr, 0 <C 1. Note that in the limit a ^ 0 and 0 ^ 0, we obtain w{x,y) = wo{x — y), a distance- 
dependent even function. However, in the fully general case (a > 0 and/or 0 > 0), it is straightforward 
to see that the function w{x,y) need not be distance-dependent or even symmetric. In particular, when 
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Fig. 2 Velocity-integration with moving bump attractors in the neural field model (2.1) with a Heaviside firing rate function 
(2.8) with threshold ^ = 0 and a cosine base weight function wo{x) = cos{x). (A) Bump of neural activity u{x,t) perfectly 
integrates velocity inputs in the case of no heterogeneity (cr = 0) and no noise (e = 0), showing the animal’s true position 
(solid line) is perfectly tracked by the center of mass of the bump (dashed). (B) In a heterogeneous network (cr = 0.1 with 
Wu = sin (a:)), the bump initially moves too fast due to a discrete attractor of the input-free system at x = 7r/2, so the 
bump’s center of mass is mismatched with the true position of the animal. (C) In the presence of spatiotemporal noise 
(e = 0.2 with cosine correlations C{x) = cos(a:)), the bump wanders diffusively so the encoded position tends to slowly 
distance itself from the true position. Here the external velocity input is constant v{t) = 0.1. Numerical simulations are 
performed using an Euler timestep with dt = 0.1 and a trapezoidal rule on the integral with dx 0.003. 


cr > 0, discrete attractors form in the network (2.1) whereby bumps tend to drift away from their initial 
position to a finite number of linearly stable locations (Itskov et al, 2011; Kilpatrick and Ermentrout, 
2013; Zhang, 1996). We consider this to be a major source of error in the network, since near-perfect 
integration of the velocity inputs could be achieved if w{x,y) = wo{x — y). For ease of analysis, the 
translationally symmetric function is typically taken to be a cosine wo{x) := cosx. We will allow Wu to 
be arbitrary, by representing it as a series of N Fourier modes 

N 

Wu{x) ;= y; an cos{nx) + sin(na;); («„) = (/3„) = 0; {al) = {^1) = al- (2.3) 

n=l 

The coefficients an^Pn are random variables drawn from the normal distribution with mean zero and 
variance cr^. 

Velocity inputs are represented by the shifting function Wy{x — y) := —Wq{x — y) as in the original 
head direction system model (Zhang, 1996) and recent grid cell models (Burak and Fiete, 2009). In the 
absence of any heterogeneity or asymmetry, the sum w{x,y) v{t) • Wy{x — y) would be translation 
symmetric but not even symmetric in general. This asymmetry produces a moving bump as the solution 
to (2.1) that will move at a speed given by \v{t)\ (Fig 2A). Incorporating heterogeneity, cr > 0, the 
system is no longer translation symmetric, and a moving bump will not move at the same speed as 
the velocity input \v{t)\ (Fig. 2B). Thus, assuming a sensory mechanism for correcting the place cell’s 
encoded position when a cue is encountered, take the velocity input to be 

v{t) := v{t)+Vc{t), 

the sum of the animal’s true velocity v{t) and an external control signal Vc{t). This is meant to account 
for the improved place representation observed when animals can employ information about sensory 
landmarks (Aikath et al, 2014; Battaglia et al, 2004; Ulanovsky and Moss, 2011). As shown in the 
schematic in Fig. lA, we assume there is a network that can access the place cell network’s perceived 
position A{t) via a readout of the center of mass of neural activity (Deneve et al, 1999) 

^(^) = [ xf{u{x,t))dx. (2.4) 

J —TV 

The present positional error is then computed by comparing the perceived position A{t) to the animal’s 
actual position given by a time integral of the velocity input 

^T{t) = [ i;(s)d5, 

Jo 
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so the error 


r{t) = Arit) - A{t), (2.5) 

which will be positive (negative) if the estimated position is to the left (right) of the true position. Note, 
we extend the domain x G [—tt, tt] to compute (2.5) in cases where the closest distance between At and 
A is across the boundary cuts at x = ± 7 r. The error r(t) is then translated either into a continuous 
velocity control signal 


Vc{t) = Ar(t) = A • {AT{t) - A{t)), (2.6) 

or a discrete control signal given by 

(2.7) 

k=l 

where sensory cues occur at times tk and A and r determine the strength and time decay of control. As 
we will show, in the case of continuous control (2.6) strengthening the sensory feedback A always leads to 
a reduction of the error. This is not the case for discrete control (2.7), since the previous sensory cue at 
tk < t becomes less relevant as t increases toward One of the main goals of this study is to explore 
how the spacing between subsequent cues — tk determines how strong A the control signal should 
be. 

The nonlinearity / is a firing rate function taken to be sigmoidal (Wilson and Cowan, 1973) 

/('^) 

where 7 is the gain and 0 is the firing threshold. For ease in analysis, we will often consider the high gain 
limit 7 ^ 00 so that / becomes a Heaviside step function of the form 


f{u) := H{u - k) 


1 : u>0, 

0 : otherwise. 


( 2 . 8 ) 


Lastly, we also will consider the impact of the additive noise increment dW{x,t). Spatially extended 
Langevin equations of the form (2.1) have become a common model of the effects of fluctuations in large- 
scale neuronal networks (Bressloff, 2012 ). The noise term is a spatially filtered spatiotemporal white 
noise process 

dW (x, t) := f T{x- y) dY{y, t)dy, 

J —TT 


where T is the spatial filter and dY{x,t) is a spatially and temporally white noise increment. With 
these definitions, the mean and variance can be calculated as {dW{x,t)) = 0 and {dW{x,t) dW{y,s)) = 
C{x — y)S{t — s), where S{t) is the delta function and C is the spatial correlation given in terms of the 
filter as ^ 

C{x -y)= f T{x- x')T{y - x')dx'. 

J —TV 

As an example, consider T{x) = cos(x) + sin(x). Then, the spatial correlation C can be computed 
explicitly as 


/ TT PTV 

T{x — x')T{y — x')dx' = / (cos(x — x') + sin(x — x')) (cos(^ — x') + sin(^ — x'))dx' 

-TV J —TV 

= / cos^(x') cos(x) cos(^) + sin^(x') sin(x) ^hi{y)dx' = 7r( cos(x) cos(^) + sin(x) sin(^)) 

J —TV 

= 7rcos(x — y) =: C{x — y). 


As we demonstrate, the control introduced to account for the impact of synaptic spatial heterogeneity 
can also be utilized to decrease errors brought about by spatiotemporal noise (Fig. 2C). 
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Fig. 3 Discrete control reduces error in imperfect velocity-integrating networks. (A) Activity plot u{x^ t) shows perfect 
velocity-integrating network is unaffected by discrete control given by (2.7), since it accumulates no errors. (B) Spatially 
heterogeneous network with cr = 0.1 and Wu{x) = sin(a:) is corrected by discrete control so that the represented position 
(dashed line) is much closer to the true position (solid line) than in the uncorrected case (compare with Fig. 2B). (C) 
Network perturbed by noise (e = 0.2 and C{x) = cos(x)) has its represented position corrected by discrete control (compare 
with Fig. 2C). Discrete control (2.9) is given at intervals tfc+i — = At = 2 with strength A = 1 and decay timescale 

T = 1. Other parameters and numerical simulations are as in Fig. 2. 


We demonstrate the impact of discrete control (2.7) on the dynamics of neural fields that imperfectly 
integrate their velocity inputs, comparing to a perfectly integrating network for reference (Fig. 3A). 
Integrating (2.7), we find that the discrete perturbations to the velocity signal are given by a series of 
exponentially decaying impulses 


Nc 

'Vc{t) = - tk). (2.9) 

fe=i 

Furthermore, note in the limit of the timestep — tk between subsequent cues tending to zero and 
r ^ 0, the equation for discrete control (2.7) approaches that of the continuous control (2.6). Assuming 
for demonstration that cues are spaced in such a way that an animal encounters one every 2 units of 
time (t/c+i — tk = 2 for all k = 1,..., Nc — 1), we demonstrate in a single realization that a network with 
spatially heterogeneous coupling can recover its ability to correctly integrate velocity (Fig. 3B). In a 
similar way, networks with additive noise can have their velocity integration corrected by the sensory 
feedback signal given by discrete control impulses (2.9) as shown in Fig. 3C. In the next section, we 
will analyze the impact of heterogeneity and noise on the position of the bump in a low-dimensional 
approximation of the bump’s center of mass A{t). 


3 Analysis and low-dimensional reduction of bump solutions 

To develop an understanding of the impact the sensory feedback signal has on the statistics of bump 
position in (2.1), we derive a low-dimensional approximation of the model that projects the dynamics 
down to a single equation describing bump position A{t). Our analysis is adapted from recent studies 
of stochastic neural field equations assuming the impact of perturbations to a translationally symmetric 
neural field can be separated into slow timescale changes to the position of patterns along with fast 
timescale changes to the profile of patterns (Bressloff and Kilpatrick, 2015; Bressloff and Webber, 2012; 
Kilpatrick and Ermentrout, 2013). Such analysis must begin by constructing the patterns of the unper¬ 
turbed system. In our case we take the velocity inputs, heterogeneity, asymmetry, and noise all to be 
perturbations to a translationally symmetric system (taking v = a = 0 = e = 0 in (2.1)). While it is 
possible to develop exact analytical results in the case wherein we break the symmetry of this model, 
which we show, it is also helpful to collect the effects of all the possible perturbations to (2.1) into a 
single scalar equation. Doing so makes it more clear how such perturbations interact. 


3.1 Stationary bump solutions to the translation symmetric network 

We begin by assuming the homogeneous connectivity function wo{x) in (2.2) satisfies evenness {wo{x) = 
Wo{—x)) and there is no heterogeneity (cr = 0) or asymmetry (0 = 0) in w{x,y). In this case, it is 
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possible to show there is a stationary bump solution u{x,t) = U{x) with U{x) > 0 over an excited 
region x G [ai,a 2 ] (Amari, 1977; Ermentrout, 1998) in the absence of velocity inputs {v = 0) in ( 2 . 1 ). 
Furthermore, the weight function is translation symmetric since 


wo{{x s) - {y + s)) = wo{x - y ^ S - s) = wo{x - y), (3.1) 

so there will be a continuum of bump locations associated with any single bump solution to (2.1) in this 
case. Stationary bumps satisfy the equation 

U{x) = r wo{x - y)f{U{y))dy. (3.2) 

J —7T 

Note that U{x + s) will also be a solution for any s, since 

U{x + s)=f woix-y)f{U{y + s))dy. 

J —7T 

A change of variables y ^ s z yields 

U{x-\-s)= [ wo{x s — z)f{U{z))dz, 

J —7T 

and another change of variables x + s x' yields 

U(x') = r wo{x' - z)f{U{z))dz, 

J —7T 


which is precisely the equation (3.2). Now, taking the high gain limit 7 ^ 00 , we employ the Heaviside 
firing rate function (2.8). Doing so allows us to generate an equation for the bump width d = a 2 — ai as 
in Amari (1977). In this case, the equation (3.2) becomes 

rCL2 

U{x)= wo{x-y)dy. 

J ai 

We then use the threshold crossing conditions U{ai) = U{a 2 ) = 0 and evenness of wq{x) to derive 


na2 na 2 —CLi nd 

U{ai)= / wo{ai-y)dy= / wo{-z)dz = / wo{z)dz =9 

J a\ JO JO 

pa2 rO pd 

U{a 2 )= / wo{a 2 -y)dy=- wo{z)dz = / wo{z)dz =9. 

J a\ J a2—ai JO 


Note that the evenness of wo{x) allows us to manipulate the above equations so they are the same 
equalities. If evenness did not hold, the above pair of equations would each be different and we would 
have an overdetermined system for the bumpwidth d, meaning stationary bumps do not exist. Thus, 

pd 

W{d)= wo{x)dx = 9 ^ d = W-^{9). 

JO 

For example, in the case of a cosine weight wo{x) = cos(x), we have 

pd 

W{d)= / cos{x)dx = sm{d) = 0 ^ d = sin ^d, 7 r —sin ^ 0. (3.3) 

Jo 

As mentioned, the threshold conditions specify the width d of the bump. Translation symmetry allows 
the position of the bump to be anywhere x G [—7r,7r], which allows this network to integrate and store 
velocity inputs as a position memory. As mentioned in section 2 , the position of the bump will be given 
by its center of mass (2.4), which for unimodal and even symmetric bumps will also be given by the peak 


A = argmax^l/(x). 


(3.4) 
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For example, in the case of cosine weight functions wo{x) = cos(x), there is an even symmetric solution 
such that ai = —a and a 2 = a, so 

r 

U{x) = / cos{x — y)dy = 2 sm{a) cos{x). (3-5) 

J—a 

Thus, the location of the bump as computed by (3.4) is Z\ = 0. Similarly, if we compute the center of 
mass using (2.4), we find 


Z\ = 


f 

J —TT 


xf(U{x))dx= / xdx = {)^ 


f" 


which is consistent. 


3.2 Perfect velocity integration by traveling bumps 

Now we explore the impact of velocity inputs {v{t) ^ 0) on the translationally symmetric network 
(w{x^y) = wq{x — y)). For now, we assume constant velocity inputs, v{t) = vq. Assuming the bump 
subsequently moves at a constant speed c, we look for a traveling wave solution u{x,t) = U{^) where 
^ = X — ct. We will show that the traveling wave speed c is exactly equal to the velocity input amplitude 
vq, under the assumption that Wy{x) = —Wq{x) in (2.1). Plugging these conditions into (2.1), we find 

-cC/'(C) + f/(C) = [ [wo{^-y)+voWy{^-y)]f{U{y))dy. 

J —TT 

Now plugging in our requirement that the velocity portion of the weight function Wy{x) = —Wq{x)^ we 
have 


-cU'iO + U{0 = -VO r - y)f{U{y))dy + T - y)f{U{y))dy. (3.6) 

J—TT J—TT 

Under the assumption that the function l/(^) satisfies the equality (3.2), we have differentiate this 
equation to yield 


V(e = r - y)f{U{y))dy. (3.7) 

J —TT 

Canceling the (3.2) portion of (3.6), we find that 

cU'iO = VO r w'oii - y)f{U{y))dy. (3.8) 

J —TT 

The equality (3.8) follows from (3.7) as long as we set c = vq. Another implication of our analysis is that 
the shape of the bump U{^) will be the same no matter what c (equivalently i;o) is, suggesting there will 
be no relaxation time if the external drive vq were to be changed abruptly. In this way, we can expect 
the translation symmetric version of the network (2.1) to integrate inputs perfectly as was originally 
proposed by Zhang (1996). 


3.3 Imperfect integration due to heterogeneity, asymmetry, and noise 

Now that we have explored the dynamics of the perfect velocity-integrating network, we study the im¬ 
pact of introducing heterogeneities (a), asymmetry (0), and noise (e) into the network (2.1). Rather than 
deriving exact solutions as we did for the translationally symmetric system, we take a perturbative ap¬ 
proach under the assumption that alterations to the symmetric system are weak. Following perturbation 
methods originally developed for the study of front propagation in reaction-diffusion systems (Panja, 
2004; Sagues et al, 2007), we employ a separation of time scales to decompose these effects into a slowly 
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evolving displacement A{t) of the bump from its uniformly translating position and perturbations to the 
bump profile This yields the following decomposition 

u{x^ t) = U{x — A(t)) + e^{x — A(t)^t) + (3.9) 


where we assume a, 0, h(t) ^ 0{e). Plugging the ansatz (3.9) into (2.1) and expanding in powers of e, 
we find that at 0(1), the system has a stationary bump solution given by (3.2). At linear order 0(e), we 
find the following equation 


ed^{x,t) =ejC^{x,t)dt ^ U\x)dA{t) ^ a f Wu{y A{t))wQ{x — y)f{U{y))dydt (3.10) 

J —TV 

-{v{t) + (p) f WQ{x-y)f{U{y))dydt + edW{x,t) 

J —7T 


where £ is a linear functional given by 


Cp{x) 



y)f'(U{y))p{y)dy 


and its adjoint operator 


*q{x) =-q{x) + f {U{x)) f wo{x - y)q{y)dy. 

J —7T 


To ensure a solution to (3.10), we require that the inhomogeneous portion of the equation be orthogonal to 
the nullspace of the adjoint operator £*. Indeed, the nullspace of £* is spanned by (fi{x) = f'(U{x))U'{x), 
where U{x) is defined by (3.2), which we can verify using integration of parts 

C*(f{x) =-(p{x) + f{U{x)) [ wo{x - y)ip{y)dy 

J —TV 

= f{U{x)) ( - U'{x) + r wo{x - y)f{U{y))U\y)dy) 

= f'{U{x)) ( - U'{x) + f_ ^ {Mx - y))fiU{y))dy) 

= f{U{x))(^- U'{x) + woix-y)f{U{y))dy'j'j = 0. 

The last line holds by differentiating the bump existence equation as in (3.7). Now, by taking inner 
products of the null vector p with the 0{e) equation (3.10), we can derive an evolution equation for 
Z\(t), the position of the bump 

/ TV nTV nTV 

(p{x)U'{x)dxdA{t)=a ^pix) Wu{y + A{t))wo{x - y)f{U{y))dydxdt (3.12) 

-TV J —TV J —TV 

/ TV pTV 

Lp{x)U'{x)dxdt ^ e / ip{x)dW{x,t)dx, 

-TV J —TV 

where we have applied the equation (3.7). We can simplify the equation (3.12) further by isolating dA(t) 
to yield the stochastic differential equation 


dA{t) = [F(Z\(t)) + v{t) + Vc{t) + 0] dt + dW{t), 
where the impact of synaptic spatial heterogeneities is described by the nonlinear function 

/y f\U{x))U'{x) fZ^Wu{y + A)wo{x - y)f{U{y))dydx 


F{A) = -a- 


r-J'{U{x))U'{x)Hx 


(3.13) 


(3.14) 
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Fig. 4 The low-dimensional equation (3.13) for the bump position A(t) provides an accurate approximation of the bump’s 
movement in the full neural field model (2.1). (A) Perfect integration of the constant velocity input v{t) = 0.05 leads to 
a constantly drifting bump (solid line) whose position is well approximated by the projected variable A(t) (dashed line). 
Inset shows the tilted potential V(A) resulting from the constant velocity input. Circles provide corresponding locations 
between the two plots at t = 10,20,30. (B) Spatial heterogeneity Wu{x) = sin(a:) with a = 0.1 causes bumps to drift 
toward local attractors of the network. Inset shows potential with a local minimum to which the trajectory is attracted. 
(C) Spatial heterogeneity Wu{x) = sin(6a:) with a = 0.2 leads to a more rapid oscillation in the trajectory A(t). (D) Spatial 
heterogeneity Wu{x) = sin(4a:) + cos(8a:) with cr = 0.1 leads to a less regular deviation in the trajectory A{t). Heaviside 
firing rate function (2.8) has threshold 0 = 0. Numerical simulations are run using the same parameters as in Fig. 2. 


and the noise term has been projected to a temporal white noise process W(t) with mean zero ((W(t)) = 
0) and variance (W(t)^) = Dt with associated diffusion coefficient 

^ 2 1-. I-. f{U{x))U'ix)f{U{y)W{y)Cix - y)dydx 

D = e - 2 -• (3.15) 

f:^f(U(x))U'(xrdx 

Setting Vc{t) = 0 and v{t) = vq (constant), the dynamics of the position variable A{t) can be equivalently 
described by a potential function 

V(Z\) = - J [FiA) +vo + (l)]dA = - j F{A)dA - {vq + (3.16) 

SO A{t) will descend the gradient of V{A) toward its local minima. Note that in the case F{A) = 0 = 
W{t) = 0, the control term will vanish Vc{t) = 0 and the bump will perfectly integrate the velocity 
input, A{t) = Jq v{s)ds. We find that the low dimensional approximation is in excellent agreement with 
simulations of the full system in this case of perfect integration (Fig. 4A). 

Ignoring the control for the time being in (3.13), we can also identify how different network imper¬ 
fections contribute to the resultant error in path integration. To do so, we simply compute the error 
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function r{t) = Z\t(^) ~ ^{t) as given in (2.5). First, note that in a network with asymmetry 0 7 ^ 0 and 
no heterogeneity F{A) = 0, the long term error accumulates linearly in time 


v{s)ds — / [v{s) + 0] ds = — 

Jo 

so the animal’s true position At will be behind (in front of) the estimated position A when 0 > 0 (0 < 0). 
We will demonstrate the impact external control via sensory cues has upon this error in section 4. Errors 
due to arbitrary heterogeneities are more difficult to express explicitly. In general, we can express the 
solution to (3.13) implicitly in this case if we assume the velocity is constant v{t) = vq: 

If indeed the function G{A) is invertible, then we can express A{t) = G“^(t), so 

r(t) = AT{t) — A(t) = V{)t — (3.18) 


r{t) = Arit) - A{t) = 


We will demonstrate some cases where we can perform this calculation explicitly in subsection 3.4. Note 
that the main impact of heterogeneities is to establish a finite number of discrete attractors, in the 
velocity input-free system, so that bumps drift toward these locations (Fig. 4B). Even in the velocity- 
driven network, spatial heterogeneities lead to a sinuous trajectory of the bump that is mismatched to a 
straight integration of velocity input (Fig. 4C,D). Lastly, note that the impact of noise can be quantified 
by averaging across realizations of the stochastic process 

A{t) = f v{s)ds-\-yV{t) = Arit) 

Jo 

While the mean position will be the same for the true and encoded positions {AT{t) — {A{t)) = 0), the 
variance will grow linearly in time 

{r{tf) = {{ATt) - A{t) f) = mtf) = Dt, 

where the diffusion coefficient D can be computed using (3.15). Previous work has characterized the 
impact of the bump profile and spatiotemporal noise correlation structure on the diffusion coefficient 
D, providing some explicit calculations (Kilpatrick and Ermentrout, 2013). In general, the main effects 
of noise perturbations on the bump will be experienced by the bump edges, where the activity variable 
u{x,t) crosses the firing rate threshold 0. We now provide some explicit calculations demonstrating the 
impact of spatial heterogeneity on the long term position of the bump. 


3.4 Explicit results for spatially heterogeneous networks with a Heaviside firing rate 

Several previous studies have characterized the impact of periodic microstructure on the propagation 
of waves in neural media (Bressloff, 2001; Coombes and Laing, 2011; Kilpatrick et al, 2008). Typically, 
periodic heterogeneities can slow down waves and even cause propagation failure. We extend these 
previous results here, showing that the low dimensional equation (3.13) allows us to estimate the location 
of bifurcations separating detectable and undetectable constant velocity inputs v{t) = vq. Again, we are 
ignoring the impact of control at this point, studying its effects in more detail in section 4. To allow for 
fully general weight heterogeneities, we consider the decomposition given by (2.3). Thus, we can integrate 
each of the Fourier modes independently to derive the function F{A) given by (3.14). Furthermore, we 
assume a cosine for the homogeneous weight function wo{x) = cos(x) and a Heaviside firing rate (2.8). 

To begin, note that the bump solution is given by (3.5) and the half-width is a = ^ [tt — sin“^ 0] as 
given by (3.3). Therefore, the spatial derivative U'{x) = —2 sin(a) sin(x). Furthermore, the null vector 
defined by (3.11) is spanned by the difference of delta distributions 6{x a) — 6{x — a). This means that 
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the frequency n cosine Fourier components of the heterogeneity, with scaling contribute the function 
F{A) in the following way 


^ {6{x + a) — 6{x — a)) cos(n(y + A)) cos(x — y)dydx 
2 sin a {5{x + a) — 5{x — a)) ^hixdx 
+ a) — 6{x — a)) (ncos(a) sin(na) — sin(a) cos(na)) cos{x)dx cos{nA) 

2(n^ — 1) sin^ a 

+ n) — S(x — a)) (cos(a) sin(na) — nsin(a) cos(na)) sin(x)(ix sin(nZ\) 

2(n^ — 1) sin^ a 

ncos(na) — cot(a) sm(na) , , 

= - ^21 -^ sin . 

In a similar way, we can compute the coefficients arising from the sine Fourier components with scaling 
(dn ^S 




FpX^) = 


{d{x + a) — S{x — a)) sm(n(y + A)) cos(x — y)dydx 
2 sin a {5{x + a) — 5{x — a)) sinxdx 
{d{x + a) — 6{x — a)) (ncos(a) sin(na) — sin(a) cos(na)) cos{x)dx sm{nA) 

2(n? — 1 ) sin^ a 

J1 + a) — 6{x — a)) (cos(a) sin(na) — n sin(a) cos(na)) sm{x)dx cos(nZ\) 

2(n^ — 1 ) sin^ a 


cot(a) sin(na) — n cos(na) , 

=-^ cos nZ\ . 

— 1 

Thus, we can write the resultant heterogeneity in general as 

N 

F{A) = a'^Cn [an sin(nZ\) - (in cos(nZ\)], 

n=l 


(3.19) 


where 


ncos(na) — cot(a) sin(na) 
— 1 


(3.20) 


and notice in the special case n = 1, we can take the limit n ^ 1 to find 

^ sin(a) cos(a) — a 
^ 2 sin(a) 


(3.21) 


We can explicitly compute the solution to (3.13) in some special cases of the heterogeneity T(Z\), 
defined by (3.19). In particular, we focus on a single cosine-shaped heterogeneity so that am = 1, o^n = 0 
(n 7 ^ m), and = 0 (Vn). Furthermore, we assume a constant input velocity v{t) = uq, so that the 
scalar equation (3.13) for A{t) is given by 

A{t) = Ksm{mA) + uq, (3.22) 


where k = aCn- Now, assuming /\(0) = 0 the equation (3.22) can be integrated to yield an explicit 
solution 


2 

A(t) = — tan“^ 
m 


yJvQ — tF tan j 

^tan ^ 

K 


V^o - 

'^0 


(3.23) 


There is a partition in (k^vq) parameter space given by the equation |/^| = |uo|. When \k\ > |uo|, 
so the arguments k,‘^ — of the square roots in (3.23) are positive, then there is a family of fixed 
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Fig. 5 Spatial heterogeneity slows and even stops the propagation of velocity-driven bumps. (A) Bump position A(t) 
demonstrates the variant propagation velocity occurring when heterogeneity (cr = 0.5) is introduced as opposed to the 
homogeneous network (cr = 0). For strong enough heterogeneity (cr = 1), propagation fails. Theory (solid line) given by 
(3.23) is well matched to simulations (dashed line) of the full model (2.1). (B) Bump velocity vi, decreases as a function 
of heterogeneity strength cr until propagation failure occurs at cr = (m^ — l)|no|/| cot (a) sin (ma) — mcos(ma)|. Firing rate 
function is Heaviside (2.8), heterogeneity is Wu(x) = cos(mx), and other parameters 0 = 0.5 and m = A. 


points of the equation (3.22), so that the bump position A{t) will eventually become pinned to a single 
position. In this case, velocity inputs are not detectable by the network, since they do not result in 
the propagation of a bump. Consistent with this, the equation (3.23) has a defined limit at t ^ oc. 


A{t) 


— tan 


— k/vq . The general formulas for all equilibria of (3.22) are given by 


2k'K 1 . A VQ 

Ak+ = -^-sm —, 

mm K 


k = 


Ah- = 


(2k + l)7r 


m 


1 . -1 1^0 
— sm —, 
m K 


0,..., m — 1, 
k = 0 ,..., m — 1. 


(3.24) 

(3.25) 


On the other hand, when |/^| < |i;o |5 the heterogeneity F{A) will not lead to pinning of bumps, so bump 
will propagate indefinitely in response to velocity inputs. However, the heterogeneity will ultimately 
reduce the speed of propagation of bumps, as found in previous studies of periodically heterogeneous 
neural fields (Bressloff, 2001; Coombes and Laing, 2011; Kilpatrick et al, 2008). We can determine the 
average reduction in the bump’s speed by calculating the time T at which A(t) crosses A = 211 jm^ 
completing one period of the heterogeneity sin(mZ\): 


T = 


27r 




which means that the average speed of the bump v\) is given 

271 jm 


Vh = _ = \lvl- K? 


(3.26) 


(3.27) 


similar to the speed scaling formulas found in Bressloff (2001); Coombes and Laing (2011). This allows 
us to directly compute the curve in parameter space at which wave propagation failure occurs, |'r’o| = 
as stated above. Note that the bump speed (3.27) depends on the frequency and amplitude of the 
heterogeneity through the term 


hv = a 


cot(a) sin(ma) — mcos(ma) 
— 1 


Thus, we can approximate coarse time-average error in path integration for this network as 

r{t) = Arit) - A{t) fa (vq - 


-K^ \t. 


(3.28) 


(3.29) 
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A more precise estimate is obtained by using the formula for A{t) given by (3.23). We demonstrate the 
accuracy of this full approximation in Fig. 4C,D. Bump position approximations (3.23) and average 
speed approximation computed from (3.27) are compared with the full neural field model (2.1) are given 
in Fig. 5. We find the low-dimensional approximation (3.22) is in excellent agreement with simulations. 
In particular, the points in parameter space at which propagation failure occur are well matched, and 
the sinuous trajectory of the bump is well tracked by our low dimensional theory. This suggests we can 
gain many insights concerning the full model by analyzing this simpler equation (3.13). 

As noted above, the existence of spatial heterogeneities in networks can lead to pinning or a reduction 
in speed of propagating bumps, which should be accurately tracking velocity-input. However, several 
previous experiments have suggested that sensory feedback is incorporated into the neural representation 
of spatial navigation (Battaglia et al, 2004; Hardcastle et al, 2015; Ulanovsky and Moss, 2011; Zhang et al, 
2014). As discussed in section 2, we propose a simple external control mechanism that incorporates a 
comparison of an animal’s current estimate of position with an external sensory cue (Fig. lA). In section 
4, we will demonstrate the improvement in position encoding afforded by sensory feedback control. 
Furthermore, we will show that there is an optimal weighting and timescale of control feedback when 
sensory cues appear discretely in space or time. 


4 Incorporating sensory cues with online control 

Recent experimental studies have shown that the presence of sensory landmarks reduces the size of 
mammalian place fields as compared to the case of no landmarks (Aikath et al, 2014; Battaglia et al, 
2004; Zhang et al, 2014). Interestingly, such a reduction in place field size can occur quite quickly, in 
response to the temporary presence of sensory information, as show in echolocating bats (Ulanovsky 
and Moss, 2011). Note here, we are referring to sensory information beyond the animal’s proprioceptive 
experience of its own motion. Specifically, we are referring to objects placed along the track of navigation 
that may provide visual, auditory, or olfactory feedback (Fig. IB). This suggests an online interaction 
between the sensory system and the path integration system that passes some positional information 
acquired by sensory cues (Tsao et al, 2013). We suggest that such a mechanisms could counteract errors 
in position encoding that could arise due to synaptic heterogeneity (Brody et al, 2003; Hansel and Mato, 
2013; Itskov et al, 2011) or noise (Burak and Fiete, 2012; Compte et al, 2000; Laing and Chow, 2001). 
However, when cues occur discretely in space, tuning the strength of feedback introduces a tradeoff 
between the immediate benefits of recent cues and the deleterious influence of older irrelevant cues. We 
explore this in the low-dimensional model (3.13) derived in section 3. 


4.1 Reducing error due to network asymmetry and heterogeneity 

We first examine the case of instantaneous cues and updates, modeled as a continuous update to the 
position, as described by (2.6). This would be the case in which landmark cues are continuously apparent 
to an animal, allowing the sensory system to send a constant stream of information to the position 
encoding network. For the time being, we also ignore the impact of noise, exploring its effect in subsection 
4.2. Under these assumption, the low-dimensional equation for bump position is 

A(t) = T(Z\(t)) (p v(t) + A(Z\ 7 ^(t) — Z\(t)). (4-1) 

As a simple example of the impact of the control term in (4.1), we examine the case of no heterogeneity 
F(Z\) = 0 and non-zero asymmetry 0 > 0. In this case, we can analytically calculate the long term 
trajectory of the error r(t) = Arit) — A(t). To do so, we can write down the first order differential 
equation for the error (Slotine and Li, 1991) 

r{t) + \r{t) = AT{t) — A{t) -|- X^Arit) — A(t)) 

= v{t) - (p- v{t) - X{AT{t) - A{t)) -^r X{AT{t) - A{t)) 

= -</>. (4.2) 

It is straightforward to calculate the solution to the linear differential equation (4.2) in the case r(0) = 0, 
finding r(t) = —p/X. Thus, perfect convergence of the trajectory A{t) to AT{t) can only be obtained 
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in the limit of infinitely strong control A ^ oo. It is also important to note that as long as the control 
strength is positive A > 0, the error r{t) will be bounded in the long time limit t ^ oo. 

We can extend our analysis of the equation (4.1) to the case of arbitrary heterogeneities using regular 
perturbation theory. Writing the linear expansion of A{t) = Z\o(t) + Z\i(t)/A under the assumption 
A ^ 1, we find that the zeroth order equation for Ao{t) is simply given by Ao{t) = Arit). Extending to 
the first order equation in 1/A, we find 

Arit) = F{AT{t)) + (p + v{t) - Ai{t). 

Applying the equation AT{t) = v{t), we thus find that Ai{t) = F{AT{t)) + 4>, which means that the 
long term error can be approximated by 


r(t) = vis)ds^ + </. + ©(l/A^) 

to first order in 1/A. Thus, as long as F{A) is a bounded function, then the error will remain bounded, 
reaching a maximum amplitude of maXa,|F(x) + 0| (Slotine and Li, 1991). 

Thus far, we have considered the case of a continuous flow of sensory information providing an 
accurate estimate of an animal’s position in space. However, in more realistic scenarios, animals receive 
external sensory information discretely in time via local landmarks (Battaglia et al, 2004; Tsao et al, 
2013) or echolocation (Ulanovsky and Moss, 2011). Sensory cues that provide a landmark for an animal’s 
present position may be captured periodically in time or more randomly; we account for both forms of 
sensory cue acquisition. As discussed in our formulation of the model in section 2, we assume the influence 
of sensory cues weakens as time elapses from the time at which they were received. This is consistent 
with recent observations concerning the evolution of place fields in bats as a function of the time since 
the last echo signal (Ulanovsky and Moss, 2011). Thus we consider the following model combing path 
integration with sensory cues acquired at times 

A{t) = F{A{t)) + 0 + v{t) + 'r’c(t), 

Nc 

Veit) = -Vc{t)/T + X^r{tk)5{t - tk), 
k=l 

= ^T{tk) — 4\(t/e). 

Analogous to the continuous control case, the error term computes the instantaneous difference be¬ 
tween the true position Arit) and the encoded position A{t) at time tk- This is then incorporated into 
the discretely incremented control term Vc{t) with strength A, and the temporal decay of cue influence is 
determined by the timescale r. We will demonstrate that for any given r, there is an optimal strength of 
feedback that trades off the error reduction of recent cues {tk) with the error increase potentially arising 
for older cues (ti, ...,t/c_i). Assuming 'r’c(O) = 0 and treating the pointwise values of r{t) as constant, we 
can integrate (4.4) to yield the piecewise smooth function 

Nc 

Vc{t) = A y; - tk), 

k=l 

as we did in section 2 for the full neural field model in (2.9). Thus, adjustments in velocity are discretely 
incremented and then decay over time. Also, note in the limit r ^ 0 and ^ 0, we obtain the 

continuous control function Vc{t) = Xr{t). This can be seen by performing this limit on (4.4) and then 
integrating. 

To demonstrate the impact of discrete control in more detail, we begin by studying the case of a 
network subject only to asymmetry (F(A) = 0 and 0 > 0). Furthermore, we focus on the case of 
constant velocity input v(t) = vq, so we can write the discretely controlled position equation (4.3) as 

Nc 

A(t) = Ip + Vo + [votk - A{tk)] - 4). (4.5) 

/c=l 


(4.3) 

(4.4) 
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We can solve the piecewise smooth differential equation (4.5) recursively, integrating with a new initial 
condition A{tk) at each cue time In the initial time domain [0,ti), /\(0) = 0 and A{t) = vq( j)^ so 
A(t) = (i;o + and A{ti) = (i;o + 0)ti. On the subsequent time domain [^1,^2), we have 

A{t) = Vo 

so 

Z\(t) = (i;o + + Arr(ti) 1 — . 

In a similar way, we can solve for A{t) on [^2,^3) to find 

2 

A{t) = {vo + (p)t +XT'^r{tk) 1 - , 

fc=i 

and in general, we can express 

Nc 

A{t) = {vo + 4>}t +XT'^r{tk) 1 - H{t-tk). 

k=l 

Thus, we can express the error as a function of time 

Nc 

r{t) = -4>t - Xr'^ritk) 1 - H{t-tk). (4.6) 

fe=i 

Expressing Vk := r{tk) and focusing on the error at the cue timepoints t/^, we can write (4.6) as 

i-i 

ri =—(l)t — Xr'^^Vk 1 — . 

k=l 

Furthermore, in the case of periodically spaced cues, we can write — tk = {At), \/k, so that 

i-i 

ri =—(/) ’ I ’ At — Xr'^^Tk 1 — . (4.7) 

k=l 

Assuming that A is not too large, the discrete equation (4.7) will have a fixed point in the long time limit 
ri ^ r*, which we can compute by taking the difference between and ri and approximating ~ r*: 

z-i 

(j>-l-At- Xrr* E ^ 

k=l 

(l)-{l + l)-At- Xrr* E [l “ Q-iW-k)At/T'^ ^ 

k=l 

and we can make the approximation 0, so that ri — yields 

0 = 0 • At + Arr* ^ r"" =—cj) • At/{X t). (4.8) 

We demonstrate the accuracy of the formula in Fig. 6A,B, showing that the error remains bounded due 
to the periodic perturbations of the discrete control term. Of course, the fixed point value given by (4.8) 
is contingent on its existence and stability. In cases where either condition is violated, the error ri will 
diverge in the long time limit (Fig. 6C). Essentially, negative feedback overcorrects for the previously 
observed errors at each cue time, leading to unstable oscillations in the error. Analytically identifying the 
cases in which ri diverges would require a more thorough study of the discrete equation (4.7). Numerical 
simulations suggest there is a boundary value of A above which these unstable oscillations occur. Thus, 
the maximal value A for which the fixed point r* exists and is stable would correspond to the optimal 
control strength, all other parameters being fixed. 


n =r = - 

n+i = r = - 
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Fig. 6 Path integration error in a network with asymmetry 0 = 0.1, discretely controlled according to (4.4). (A) Error 
resulting from asymmetry plus discrete control with time decay r = 1 quickly reaches the fixed point r* (thin lines) given 
by (4.8). Notice as the strength of control A is increased, the long term error decreases. Low dimensional approximation 
(solid lines) given by (4.5) is in excellent agreement with simulations (dashed lines) of the full model (2.1). (B) Increasing 
the control decay timescale r = 2 leads to longer lasting oscillations in the error. (C) Making the control too strong, 
A = 4.5, leads to instability in the error. Negative feedback produces oscillations that grow in amplitude. Time spacing 
between cues is At = 1. Other parameters 0 = 0.5, vq = 0.1. Numerical simulations utilize the same parameters as in Fig. 
2. 


We now study the case of heterogeneities and explore the impact of sensory cues on the long term 
error. Note, in the case of no asymmetry and constant velocity input v{t) = vq, the low dimensional 
equation for bump position takes the form 


Nc 


A{t) = F{A{t)) + + A y] [vQtk - a(4)] e 




H{t - tk). 


k=l 


(4.9) 


While we cannot solve (4.9) explicitly for general heterogeneities F{A)^ we can numerically analyze 
the impact of both the control strength A and the control decay timescale r on the long term error 
r(t) = Z\t(^) — A{t). Specifically, we associate error with a scalar quantity by computing the log of the 
I/^-norm 


R := In ||4 \t(^) — ^(^)||2 = In 



I^t(^) ~ A(t)\^dt 


(4.10) 


where t/ is time at which the path ends. We compare the effects of varying the spacings tk+i—tk between 
subsequent cues, testing both time-periodic cues (t/c+i — = At, V/c) and exponentially distributed 

spacings {p{At) = Furthermore, we randomize the heterogeneity according to the formula 

(2.3) with variance = 1 with four total modes {N = 2). To average error across many realizations, 
we simulate the controlled version equation (4.9) for many different randomly generated heterogeneities, 
compute an L^-norm of error Rj for the jth realization and average (R) = ^ realizations. 

We are mainly interested in the (A,r) values that minimize the average log error (R). Our findings 
are summarized in Fig. 7. First, we note that there is always a curve through (A, r) space determining 
the values of the control term that minimize the average error (R). In all plots, the associated r value 
decreases with A and vice versa. In general, we find this relationship to be roughly inversely proportional 
A oc I/r. This means that stronger control should decay more quickly, and equivalently weaker control 
can last longer. Furthermore, by comparing plots for periodic cues with At = 4 (Fig. 7A) versus At = 2 
(Fig. 7B), we find longer decay timescales associated with each A value in the case At = 2. Such a 
trend may arise due to the fact that more frequent updates in sensory information via cues prevents 
overcorrection that could occur in the case of less frequent curves. A similar trend arises in the case of 
exponentially distributed time spacings between cues (/i = 0.5 in Fig. 7C versus /i = I in Fig. 7D). 
When cues are more frequent, the optimal timescale of decay r is slightly larger for each value of A. 

In addition, we have studied the average error as a function of time for both continuously and 
discretely controlled networks with heterogeneities. Note that we randomize the heterogeneity Wu{x) = 
ai cos(x) so that is normally distributed with variance unity. To compute the average error, we take 
the mean of the absolute value {\r{t)\) as shown in Fig. 8. In the case of strong continuous control, it is 
possible to substantially decrease the impact of heterogeneities as compared with the uncontrolled case 
(Fig. 8A). Discrete control maintains an intermediate level of error, since there is not a constant stream 
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Fig. 7 Average log error (R) computed across realizations of (4.10) for the discretely controlled low-dimensional approxi¬ 
mation (4.9) with spatial heterogeneity resulting from (2.3) with N = 2 and coefficient variance cr^ = 1. (A) For periodically 
appearing control cues with At = 4, we find there is an intermediate curve (solid) of (A, r) values that minimizes (R). In 
particular, as the control decay timescale r is increased, the optimal value of A decreases. (B) The same trend is consistent 
for periodic cues with spacing At = 2, but the curve of optimal (A, r) values shifts so there are higher values of r asso¬ 
ciated with each value of A as compared with A. (C,D) When spacings between cue times are exponentially distributed 
p{At) = with p = 0.5 in C and /r = 1 in D, we find the optimal curve shifts to shorter values of r for more frequent 

cues. Other parameters vq = 0.15, 6 = 0, cr = 0.1, and simulation time tf = 40. Numerical simulations are performed using 
Euler’s method with a timestep of 0.05, and each grid point used 1000 realizations. 


of information provided to reduce error. Varying the strength A and timescale r of control alters the long 
term variance in the error (Fig. 8B). As suggested by Fig. 7, strong and fast decaying control tends to 
lead to substantial reductions in error. 


4.2 Reducing error due to dynamic fluctuations 

We now examine the impact of sensory feedback on networks subject to temporal noise fluctuations. 
Dynamic variability in networks can arise from ion channel fluctuations (Chow and White, 1996), synaptic 
variability (Ribrault et al, 2011), or finite size effects (Bressloff, 2009). As demonstrated in our analysis 
in section 3, we can derive a reduced equation for the position of a velocity-driven bump subject to noise 
to a single stochastic differential equation (3.13). Focusing specifically on the impact of noise, taking 
constant speed v{t) = vq, and ignoring heterogeneities, we find that the controlled equation for the bump 
position takes the form 


dA{t) = [i;o + Vc{t)] dt + dW{t). 


(4.11) 


We begin by examining the case of continuous sensory feedback, in which case (4.11) becomes 


dA{t) = [i;o + Xvot — XA{t)] dt + dW{t). 


(4.12) 
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Fig. 8 Average error {|r(t)|} as a function of time in a heterogeneous network (wu{x) = ai cos{x) with control. Velocity- 
input is constant v(t) = vq = 0.05. (A) Numerical simulations (dashed lines) of the neural field (2.1) are well matched 
by the low-dimensional approximation (solid lines) given by (4.9). As demonstrated, continuous control provides the best 
reduction in error, but discrete control with At = 4 still provides an appreciable reduction. (B) Plot demonstrates the 
impact of varying the control strength A and the control decay timescale r. Other parameters are 6 = 0.2 and a = 0.1. 
Numerical simulations are run with the same parameters as in Fig. 2 for 1000 realizations each curve. 


Note that (4.12) is a non-autonomous Ornstein-Uhlenbeck process, and we can use integrating factors 
to identify an explicit solution. Utilizing the change of variables h(Z\,t) = and differentiating with 
respect to t, we find 

dh{A, t) = dAe^^ + XAe^^ = e^^ [i;o + Xvot — XA + XA] dt + e^^dW 

= e^^ [i;o + Xvot] dt + e^^dW = d (^vote^^) + e^^dW. (4-13) 


Assuming /\(0) = 0, we can integrate (4.13) and multiply through by e to yield the solution 

A{t) = vot + e-^^ [ e^^dW{s), 

Jo 

whose mean is {A{t)) = v^t and variance is given 

where the diffusion coefficient D can be calculated from the neural field model parameters as in (3.15). 

^/d 

The long term variance is thus given by lim {A{t)^) — {A(t))‘^ = ——. Note that in the limit A ^ oo, 

t^oc 2A 

the variance goes to zero {A{t)‘^) — {A{t))^ 0, suggesting that strengthening continuous control will 

always reduce the average error further. Continuous control substantially reduces the long term variance 
in the bump position A{t) as well as the error 


r{t) = Arit) - A{t) = [ 

Jo 


As 


dW{s). 


Note that (r(t)) = 0 and {r{t)‘^) = ^ [l — . We compare the continuously controlled system to 

the control-free system in Fig. 9A, revealing the long term saturation in the position variance. 

We also study the effect of discrete control on the variance in position, using the low-dimensional 
approximation of bump position 


dA{t) = [v{t) -h Vc{t)] dt + dW{t), (4.14) 

Nc 

Vc{t) = ^ X^r{tk)S{t - tk). 

k=l 
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Fig. 9 Variance {A{t)^) — {A{t))^ computed for the noise-driven network with control. Velocity-input is constant v{t) = 
vq = 0. (A) Numerical simulations (dashed lines) of the neural field model (2.1) are well matched to low-dimensional 
approximation (solid lines) given by (4.11). Continuous control substantially reduces the variance, but discrete control with 
At = 4 also provide a variance reduction. Notably, the variance saturates in the case of discrete control as well. (B) Similar 
to the case of quenched variability through heterogeneity in Fig. 8, varying the strength and timescale of control alters the 
long term variance. Other parameters are 6 = 0.2 and e = 0.1. Numerical simulations are run with the same parameters as 
in Fig. 2 with 1000 realizations each curve. 


Again, this is under the assumption that cues are spaced discretely in time or space, and they provide 
sensory input for a brief period of time. As in subsection 4.1, we can solve (4.14) iteratively. To begin, 
note that when t G [0,ti), A{t) has yet to be affected by the feedback control term in (4.14), so A{t) = 
v^t + yV{t). Subsequently, we can integrate (4.14) to find the stochastic formula for A{t) after the first 
cue at ti’. 


A{t) = vot + Xr{voti - A{ti)) 


I _ Q-(t-ti)/T 


+ W(0, 


t G [ti,t2). 


Each realization will typically produce a different value for A{ti) = voti-\-yV{ti). Note that {A{t)) = vot, 
so the impact of noise and control can be observed by calculating the variance (Gardiner, 2004) 

^ (vot + XT{voti - A{ti)) (l - + W{t)j ^ - (votf 

= {W{tf) + (^{XTfWihf (l - - (2XTW{ti)W{t) (l - 

=Dt + XrDti (l - (At (l - - 2) . (4.15) 

One insight to be gained from the formula (4.15) is that infinitely strong and fast decaying control, even 
when it is discrete, will minimize the variance in the position. Specifically, if we take A = Aq/t, then we 
can write 

{A{tf) - {A{t)f =Dt + XoDti (l - (Aq (l - - 2) . 

Taking the limit as A ^ oc, we find that 

lim [{A{t)'^) - {A{t))^] =Dt + XoDti{Xo - 2), 

which is minimized when the scaling term Aq = 1, yielding {A{tY) — {A{t))‘^ = D(t — ti). 

We can solve (4.14) explicitly for an arbitrary number of cue times, yielding 


Nc 

A{t) = V{)t + Ar E F(t/c) 

k=l 


I _ Q-it-tk)/r 




(4.16) 


While it is clear that the mean {A{t)) = vot^ it is more complicated to compute the variance (Z\(t)^) — 
(Z\(t))^ in general. This is chiefly due to the fact that r{tk) will depend on (r(ti),..., r(t/c-i)), and this 
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Fig. 10 Average log error (R) computed across realizations of (4.10) for the discretely controlled low-dimensional approx¬ 
imation (4.16) driven by additive noise with amplitude e = 0.1 as described in (2.1). (A) For periodic cues with spacing 
At = 4, we find that the curve (solid line) of optimal (A,t) values has r decreasing with A as in Fig. 7. (B) A similar 
tend is observed for exponentially distributed p{At) = spacings between cue times with p = 0.5. Other parameters 

i;o = 0.15, ^ = 0, and simulation time tf = 40. Numerical simulations are performed using the same method as in Fig. 7. 


long-lasting history-dependence will accumulate indefinitely. To gain some analytical understanding, we 
make the assumption of brief control impulses, so that r ^ 1 and 1, V/c. In this case, we 

can write the equation for the update of the error term Vk := r{tk) as 

n+i - (1 - Ar)r/e + W(4+i) - W(4), 

where ri = W(ti). Again, it should be clear there is an inverse relationship between the impact of A and 
r on the long term error in this limiting case. We compute the variance numerically from (4.16) for the 
case of discrete control in Fig. 9, demonstrating an excellent match with the neural field model (2.1). 

We conclude by computing the average log error (4.10) across realizations of (4.14) in Fig. 10. Notice 
again that the optimal value of r, which minimizes (R) is inversely related to the strength of control. Fur¬ 
thermore, this trend is preserved whether cues appear periodically in time (Fig. 10A) or at exponentially 
distributed intervals (Fig. lOB). 


5 Discussion 

We have introduced and studied a neural field model of path integration with sensory feedback. Velocity- 
input results in the propagation of a bump attractor whose position encodes an animal’s estimate of its 
position. Sensory information is assumed to come in the form of cues that are either constantly present, 
in the case of continuous feedback, or present at discrete points in time, in the case of discrete feedback. 
The full neural field model (2.1) can be reduced to a single scalar equation (3.13) for the resulting position 
of the bump attractor. Analyzing this reduced equation, we have found that continuous control can be 
used to reduce the error to zero in a variety of cases. Incorporating the more realistic assumption of 
discrete sensory control, we find a tradeoff arises as the strength of control A is tuned: error reduction 
when cues are recent, counteracted by error increases when cues are older and irrelevant. Thus, there is 
an optimal control strength A that minimizes the long term error in the model’s position estimate. This 
pattern holds when errors originate from spatial heterogeneities as well as dynamic fluctuations. 

Our analysis has focused on one-dimensional periodic systems, wherein it is assumed the animal is 
navigating along a narrow annular track (Fig. IB). This was based on the protocol used in the exper¬ 
iments of Battaglia et al (2004), which were used to study the effect of local cues on the sharpness 
of neuronal place fields. However, there are several studies of navigation in two-dimensional and even 
three-dimensional space that demonstrate mammals’ ability to use sensory cues to perform error correc¬ 
tion (Geva-Sagiv et al, 2015; Solstad et al, 2008). For instance, a recent study has demonstrated that 
encounters with the boundaries of rectangular environments correct for the systematic drift in position 
representation (Hardcastle et al, 2015). In particular, border cells in medial entorhinal cortex (MFC) 
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are thought to provide inputs to position-encoding grid cells when an animal senses an environmental 
boundary. Such recent studies are consistent with the predictions of planar models of spatial naviga¬ 
tion based on the dynamics of velocity-driven bump attractors (Burak and Fiete, 2009; Samsonovich 
and McNaughton, 1997). The model we have presented here could be extended to incorporate the ef¬ 
fects of position-dependent cues, like boundaries, in two-dimensional domains. We expect the extension 
to two-dimensional neural field models should be possible through a similar negative feedback control 
mechanism to those presented in section 2. Our derivation of the reduced equation would then simply 
yield a position variable that is two-dimensional, with a correction term along each coordinate. 

We also note that there is recent evidence that the position of discrete objects in the environment 
may be encoded by cells in the lateral entorhinal cortex (LEG) (Tsao et al, 2013). In particular, these 
cells tend to be inactive in open environments with no spatial landmarks, but they become active in the 
presence of objects that can help animals to orient themselves (Deshmukh and Knierim, 2011). Some 
cells in LEG, object-trace cells, have been shown to fire when an animal encounters a location where an 
object previous was located, demonstrating a persistent memory of location (Tsao et al, 2013). If in fact 
such cells provide inputs to the position-encoding networks in MEG or hippocampus, LEG object cells 
could provide a candidate mechanism for the sensory feedback control which we have modeled here. 
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